Magnetothermoelectric Response at a Superfluid— Mott Insulator Transition 
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We investigate the finite temperature magnetothermoelectric response in the vicinity of a 
superfiuid-Mott insulator quantum phase transition. We focus on the particle-hole symmetric 
transitions of the Bose-Hubbard model, and combine Lorentz invariance arguments with quantum 
Boltzmann calculations. By means of an epsilon expansion, we find that a non-vanishing thermo- 
electric tensor and a finite thermal conductivity are supported in this quantum critical regime. We 
comment on the singular Nernst effect in this problem. 
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Since the discovery of high temperature superconduc- 
tivity, quantum phase transitions between Mott insula- 
tors and supcrfluids/superconductors have been the focus 
of intense scrutiny. More recently, remarkable advances 
in the control of cold atomic gases have allowed the ob- 
servation of such phase transitions in systems of bosonic 
atoms . Such transitions separate two of the most fasci- 
nating highly correlated states of matter: the superfluid, 
which reveals the dramatic consequences of quantum co- 
herence on the macroscopic scale, and the Mott insulator, 
which affirms the impact of strong interactions. In par- 
ticular, a continuous quantum phase transition between 
the two clearly requires a combined treatment of strong 
quantum fluctuations and interactions. Understanding 
this interplay is useful not only to the particular systems 
at hand, but to a wide variety of condensed matter prob- 
lems 0- 

Our general objective in this present work is to ad- 
vance the theory of continuous superfluid-Mott insulator 
(S-MI) transitions in bosonic systems. In particular, we 
wish to focus on their quantum critical regimes, and to 
examine the response to applied electric and magnetic 
fields, and temperature gradients. We are motivated by 
measurements of such magnetothermoelectric response 
by Ong and co-workers 0, Q . These have yielded crucial 
insights into the importance of superconducting fluctu- 
ations in the cuprates although a direct correspon- 
dence with our present bosonic system is clearly lacking. 
Our results apply most simply to charged bosons, al- 
though there are plausible suggestions to generate pseudo 
magnetic fields in neutral atomic systems; see e.g. 0. 

Specifically, we focus on the most straightforward 
example of a S-MI quantum phase transition — the 
particle-hole symmetric version occurring at the tips of 
the Mott lobes in the ubiquitous Bose-Hubbard model 
|3| — and the quantum critical regime of this transition. 
This model arises in numerous contexts including ^He 
layers, thin film superconductors and Josephson junction 
arrays. More recently, this system has also been realized 



in ground breaking experiments on cold atomic gases , 
where it is possible to tune through the S-MI transition 
using lasers. This approach confirms the predicted phase 
diagram [3 and offers a unique handle on these strongly 
correlated systems. 

In view of this broad spectrum of applications, the 
Bose-Hubbard model has received considerable theoret- 
ical attention, see for example H S Oi Q and 
references therein. Here we shall incorporate the signifi- 
cant effect of an applied magnetic field on the quantum 
critical regime of the particle-hole symmetric transition. 
This paper is structured as follows: We begin with a short 
introduction to the Bose-Hubbard model and its quan- 
tum phase transitions. We recall the quantum Boltz- 
mann equation (QBE) and the epsilon expansion used to 
access the finite temperature critical dynamics 0|. We 
briefly discuss the response to a single external field, and 
then move on to the situation in crossed electric and mag- 
netic fields. We combine Lorentz invariance arguments 
with Boltzmann calculations to provide the fundamental 
response coefficients. We then examine the complemen- 
tary situation in a crossed thermal gradient and a mag- 
netic field. We conclude with implications for possible 
Nernst measurements, a brief summary, and directions 
for further research. 

Model — The Bose-Hubbard model describes bosons 
hopping on a lattice with amplitude t, and interacting 
via the short range repulsive interaction U: 
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The Bose creation and annihilation operators satisfy the 



b]b, 



usual commutation relations, [bi^w^j — "jjj 
the number of bosons at site z, and fj, is the chemical 
potential. In the context of a Josephson array or super- 
conductor, the bosons represent Cooper pairs of charge 
Q = 2e, tunnelling between superconducting regions. In 
general, one may also include the effects of disorder and 
long range interactions into such a model, but here we 
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shall concentrate on the simplest case Q . The phase di- 
agram of the Bose-Hubbard model exhibits both super- 
fluid and Mott insulating regions , the latter occurring 
for strong enough repulsive interactions. In particular, 
as a function of the chemical potential, yu, this model 
exhibits a series of Mott insulating "lobes" where the 
density of bosons is pinned to successive integers. At the 
tip of these lobes, the energy cost to producing particle- 
hole excitations vanishes and the model is particle-hole 
symmetric. In the vicinity of these points, the S-MI tran- 
sition is described by the continuum action of a complex 
scalar field <& 0, 



5 = 



(2) 



where D ~ d + 1, and d is the number of spatial dimen- 
sions. This is nothing but a quantum Landau-Ginzburg 
theory for the superconducting order parameter, The 
bare couplings m and uq may be related to the mi- 
croscopic parameters of the Bose-Hubbard model, and 
we refer the reader to for more details. Away from 
these particle-hole symmetric points, the effective action 
picks up an additional term linear in the time derivative. 
Correspondingly, the dynamical exponent changes from 
z=ltoz = 2|3|. For the purposes of this paper we will 
focus on the case with z = 1. 

Quantum Boltzmann Equation — In the absence of in- 
teractions, the model (O is the celebrated Klein-Gordon 
theory. The two modes of opposite charge correspond to 
density fluctuations, or particle-hole excitations, of the 
Mott insulator. Performing this mode expansion on the 
full model (0) allows one to develop a Boltzmann ap- 
proach to transport which is particularly simple at non- 
zero temperatures and small applied fields [iflj. The re- 
sulting quantum Boltzmann equation (QBE) is a nonlin- 
ear integro-differential equation for the momentum space 
distribution functions, /±(k, t), of particle and hole exci- 
tations 



dt 



± Q (E Vk X B) 



■ 9k 



I±[/+,/-], (3) 



where Vk = 9ek/9k and £k = = vk^ + Jr?. For 
simplicity we consider a spatially homogeneous system 
in uniform external fields. The collision term represents 
scattering between these excitations, and incorporates 
the nonlinear interaction of the Landau-Ginzburg action 
(0) in a systematic cpsilon expansion |lOj : 



d^'k. 



2n¥2e 



ki 



{27rf+'5ik + ki - k2 - k3)Sie + Si~S2- £3), (4) 
where scattering out of state k is given by 

.^r = 2/±(k)/^(ki)[i+/±(k2)][i+/^(k3)] 

+ /± (k) /± (ki ) [1 -f /± (k2 )] [1 + /± (k3 )] , 



and we have suppressed the explicit time dependence of 
the distribution functions. Scattering into state k follows 
by interchanging f± and 1 + /±. At the Wilson-Fisher 
fixed point, where the renormalized mass vanishes, the 
bare couplings must be tuned to the values [loj 
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uo = 



(5) 



The structure of the collision term Q may also be seen 
using Fermi's Golden rule, where the l + f factors remind 
us that we are dealing with a system of bosons. In this 
representation the electric current takes the form 



d'^k 



and the heat current is given by 
d'^k 



vk[/+(k,t)-/_(k,t)] (6) 



3h = J ^|A.Vk£k[/+(k,t) + /-(k,t)]. (7) 

In writing this last expression we have used the fact that 
within linear response, the heat current coincides with 
the energy current. The transport coefficients of interest 
are defined by the relationships 
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a a 
a K 



E 



(8) 



In order to work below d = 3 in the epsilon expansion, we 
will restrict ourselves to working with E and VT orthog- 
onal to B, and to transport coefficients in the "xy" plane 
defined by E and E x B. Particle-hole symmetry places 
constraints on these remaining coefficients: the longitu- 
dinal component of a and the transverse (Hall) compo- 
nents of a and k must vanish. Finally, we note that 
here, as in the standard Boltzmann treatment of magne- 
tothermoelectric transport, we do not explicitly consider 
magnetization currents induced by the applied magnetic 
field 0, 01 . We will discuss the implicit cancellations 
behind this elsewhere 

Separate response — Before we turn our attention to 
the general problem in combined fields, let us first discuss 
what happens when each of the fields, E, VT, and B, is 
taken separately. The application of an electric field at fi- 
nite temperatures leads to a linear response conductivity, 
axx, calculated in 0|. While the model thus supports 
current relaxation, it does not support energy current re- 
laxation; it is known that its response to an applied tem- 
perature gradient — the thermal conductivity 
finite [isj . The inclusion of the highly irrelevant umklapp 
scattering, reminiscent of the one-dimensional problems 
studied in [l^ , will ultimately render this quantity finite, 
though anomalously large as T ^ 0. Finally, let us now 
turn to the response to a magnetic field taken alone. In 
the quantum critical region, and in the absence of an 
electric field and temperature gradient, the only relevant 
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energy scales are the temperature and the magnetic field. 
On general grounds we expect the free energy density to 
scale as 



T{T, B) = T^+d/- 



B 



(9) 



where /i is some scaling function. Here we have used 
the fact that the correlation length, ^, diverges with the 
correlation time, ~ 1/?^, according to ^ ~ (6)^^^- 
The prefactor is thus an energy density. In addition we 
have used the fact that the vector potential, |A| ~ 
as may be seen from its gauge transformations, and so 
|B| - 1/2.2 ^ T^/^ It follows that the linear response 
magnetization for example, scales as 



M 



dB 



(10) 



This is consistent with a finite temperature, diagram- 
matic Kubo calculation of the magnetic susce ptib ility of 
a charged scalar field, with m ^ T and z = 1 [iJl . 

Crossed E and B fields — The motion of a single rela- 
tivistic charged particle in crossed electric and magnetic 
fields has two distinct regimes of behavior the 
regime where |E| < c|B| one may always find a frame 
moving with the velocity 



vd 



E X B 



(11) 



where the electric field vanishes. In this moving frame the 
particle executes cyclotron orbits. Boosting back to the 
lab frame the trajectories are helical, with a well defined, 
charge independent, transverse drift velocity Hll|l . In the 
regime |E| > c|B| it is possible to find a frame where 
the magnetic field vanishes, but it is no longer possible 
to make the electric field vanish. In consequence, in the 
absence of any dissipative processes, the particle is con- 
tinually accelerated by the electric field. In the lab frame 
the trajectories are hyperbolic. Since we are interested in 
linear response in a background magnetic field, wc shall 
focus primarily on the drift regime; we comment on the 
hyperbolic regime later. The speed of light at issue here 
is the bosonic mode velocity suppressed in Eqn. |(2J). 

Drift regime |E| < c|B| — So far wc have used the 
Lorentz transformations to note some properties in the 
absence of the collision term. However, the Lorentz struc- 
ture also enables us to make progress with the full QBE 
©. As discussed above, the electric field vanishes in a 
frame moving with velocity Vq. Since a magnetic field 
does not affect the energy of a particle it follows that 
a thermal distribution holds in this frame. Using the 
Lorentz transformation for energy this suggests that 



/±(k) = /o(4) = /o 



£fc - vp.k 

v/1 - «d/c2 



(12) 



is a solution of the full QBE including the collision 
term. This remarkable fact may be verified by explicit 
substitution. In particular, we may Taylor expand H12|) 
to linear order in Vq, or equivalently the electric field. 
Substituting this into iQ yields a non-vanishing trans- 
verse heat current, and a corresponding thermoelectric 
tensor 



2c' 
dBT 



d'^k 
{2nhy 



dfo 

d£k 



^ • (13) 



This result also follows from entropy drift [l^. Note 
that, we may also make contact with (|12f) from a direct 
linearization of the QBE. Substituting the decomposition 

/±(k) = /o(£fe) ± Qk.E 7/;(fc) + k.(E x B) ^^(fc) (14) 

into equation ||2Jl and linearizing in the electric field, 
yields '0(^) = 0, and an expression for ipj_{k) consistent 
with the Taylor expansion of (|12() . Having established the 
validity of H13|l as a legitimate result for the full QBE, 
we may complete the final steps of the epsilon expan- 
sion. Following Damle and Sachdev to lowest order 
we compute the numerical prefactor in d = 3 and for 
0: 



m 



_ An' kg 



he 



(15) 



The form of this result is consistent with scaling argu- 
ments in the vicinity of a QCP 



axy{T, B) 



rp{d-2)/z ^ 



B 



where /2 is another scaling function. Note that the de- 
pendence on 1/B, which stems from the drift velocity 
(tTT|l . is analogous to that of the clean, single carrier clas- 
sical Hall conductivity. It is important to recall how- 
ever, that axy = 0, in our particle-hole symmetric, two 
species case. Likewise, axx = in this regime, owing to 
the purely transverse nature of the distribution function 
((T^ . Having completed our analysis in the drift regime, 
let us discuss the transport properties in the presence of 
a thermal gradient and a magnetic field. 

Crossed VT and B fields — In keeping with our pre- 
vious discussion, let us consider |VT| <C |B|. We may 
introduce a temperature gradient into the QBE by allow- 
ing the temperature variable to be a function of position 
[iSj . The generic QBE, with E = 0, reads 



df± , df± df± 
^+v,.— ±Q(v,xB).— 



I±[/+,/-]. (17) 



In the absence of any material inhomogeneity we may 
assume that any spatial variation is due to the imposed 
temperature gradient: 



df± 
9x 



= VT 



df± 
dT 



T dek 



1& ■ (18) 
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In the last step we have used the fact that we are inter- 
ested in hnear response in VT. To solve the resulting 
linearized equation we parametrize 

/±(k) -/o(efc) + k.UV(fc)±Qk.(Ux B)V±(fc). (19) 

To lowest order in the epsilon expansion we may drop the 
collision terms to obtain 



^(fc) = o, = 



Q2|B|2 



dfo 



(20) 



This reproduces our previous result (|13|l as expected from 
Onsager reciprocity. 

Proceeding to O(e^) we find a non- vanishing longitu- 
dinal displacement 

V'(fc) = (^) &i[VJi(fc)Fi+^x(fci)F2], (21) 

where V'-l(^) is given by our zeroth order result (|20|l . 
and Fi(fc, fci) and F2(fc, fci) are the kernels employed by 
Damle and Sachdev [lOj. Introducing the dimcnsionless 
variable, k 
be written 



ck/ksT, the longitudinal displacement may 



m = 1- 



fkeT 



h \QBc 



G{k), 



(22) 



where 



G{k) = k 



dh [P(fc)$i(fc,fci) + *2(fc,fci)P(fci)] . 

_ - - (23) 
Here P(fc) = fee'"'/ {e^ — 1)^ is the massless limit of ip±{k) 
and $1,2 are the rescaled kernels Upon integra- 

tion over momentum space, the longitudinal displace- 
ment leads to a finite thermal conductivity 



(24) 



where Ib = ifi/QBY^^ and At = hc/ksT are the mag- 
netic length, and a suitable thermal wavelength respec- 
tively. To lowest order in the epsilon expansion the nu- 
merical coefficient is given by 



g = -^2 / dkk^Gik) « 5.55. 



(25) 



This equation mirrors (3.33) of Damle and Sachdev |10l| . 
The extra factor of momentum arises because we are con- 
sidering heat transport as opposed to electrical transport. 
In addition, the scaling function G{k) is determined by 
the direct integration given in equation (|23|l . It is not ob- 
tained by solving a linear integral equation such as (3.31) 

of [13. 

Beyond the linear Boltzmann regime — Above we have 
considered linear transport response in E and VT with 
a finite B <^ T already present. This is what is nat- 
urally described by the linearized Boltzmann equation. 



Evidently, this is not all the behavior that is possible in 
the quantum critical regime: with T ^ E,B,\7T^ it is 
still possible to have other ratios for the three applied 
fields. We have previously noted that for E > cB sin- 
gle particle orbits are fundamentally different. One can 
show more generally that the drift solution fails to exist 
in this domain, even in the presence of the interactions 
contained in the fixed point theory (0). It seems clear 
that the fixed point theory should yield a finite a^x , and 
an infinite axy in this regime, but a formal demonstration 
needs to be constructed. The situation with an applied 
VT is similar, but in need of more work. Specifically, the 
boundary between the VT <C B regime with a finite Kxx, 
and the VT 3> B regime which is expected to exhibit an 
infinite Kxx, although finite axy^ cannot be located with 
precision by our methods. We hope to report progress 
on these remaining issues elsewhere [l^ . 

Concluding remarks — We have examined the magne- 
tothcrnioclectric response at the S-MI transitions of the 
Bose-Hubbard model. In linear transport in an applied 
magnetic field, we obtain a finite thermoelectric tensor 
and thermal conductivity. It is worth noting that in con- 
trast to the i? = electrical conductivity [l^ these coef- 
ficients are regular in the epsilon expansion. Finally, we 
comment on a Nernst measurement in such a system in 
which we measure a transverse electric field induced by 
a thermal gradient in an open electrical circuit. For the 
measurement to be well defined, we will need to invoke 
irrelevant umklapp processes to regularize quantities that 
arc infinite in the fixed point theory. This done, we find 
that in the scaling limit, an infinitesimal VT induces a 
finite Nernst electric field of strength B/c. 
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